Constraining an imparity violating supergravity model with the Higgs induced 

Majorana neutrino magnetic moments 
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It is well known, that i?-parity violating supersymmetric models predict a non-zero magnetic 
moment for neutrinos. In this work we study the Majorana neutrino transition magnetic moments 
within an RpV modified minimal supergravity model. Specifically, we discuss the contributions 
coming from the charged Higgs bosons, higgsinos, leptons, sleptons, charged gauge bosons, and 
charginos. We use the experimental results from the MUNU collaboration to restrict the model's 
parameter space. A comparison with two other types of contributions (only trilinear RpV couplings 
and trilinear plus neutrino-neutralino mixing) is also presented. We have found that the presently 
discussed processes dominate significantly, exceeding in some cases even the experimental limits. 
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I. INTRODUCTION 



Neutrinos are one of the most mysterious particles 
known. They interact very weakly with matter, taking 
part only in the weak (and presumably gravitational) in- 
teractions. Due to mixing [l| their mass eigenstates and 
interaction eigenstates differ substantially, which gives 
rise to the generation lepton number violating oscillations 
. They are also the lightest fermions ever observed - in 
fact the smallness of neutrino masses cannot be explained 
using the standard Higgs mechanism 0. All this indi- 
cates, that neutrinos are a window to some new physics 
beyond the standard model. 

The possible electromagnetic properties of neutrinos 
in the form of an induced magnetic moment are partic- 
ularly interesting, since the electromagnetic interactions 
are much easier to control in any precise measurement. 
This problem has been studied in the standard model 
[3, Q and its extensions. In Ref. @ the magnetic mo- 
ments have been reviewed in models with spontaneously 
broken lepton number. A lot of work has been devoted 
to discuss the neutrino masses in the i?-parity violat- 
ing (RpV) supersymmetry 0-113] • I n these models, the 
neutrino mass is generated mostly at the one-loop level 
through an RpV loop containing in the simplest ver- 
sion a particle-sparticle pair (lepton-slepton or quark- 
squark). The approach usually used in the literature 
assumed no mixing among particles, and some average 
common value for the masses of the superparticles at low 
energies. In Ref. [HJ a basis independent classification 
of various diagrams leading to Majorana neutrino masses 
has been presented. The study of neutrino masses can be 
directly extended to the magnetic moment case by intro- 
ducing an interaction with an external photon. 

A more systematic approach to the calculations of the 
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magnetic moments has been presented in our previous 
works 0, [11-13. The SUGRA model has been used 
to calculate masses and magnetic moments in Ref. [llij . 
The contribution coming from neutrino-neutralino mix- 
ing has been discussed in Ref. [TtJ , while the contribution 
from generic diagrams with two mass insertions on the 
internal lines was presented in (l3| . Calculations within 
the gauge-mediated model were performed in [l8|]. We 
have studied the Majorana neutrino masses and transi- 
tion magnetic moments using exact analytic formulas and 
up-to-date experimental data concerning neutrino oscil- 
lations and the non-observability of the neutrinoless dou- 
ble beta decay (Oz/2/3). Our numerical procedure rested 
upon the assumption that most of the free parameters 
of the supersymmetric model unify in certain way at the 
GUT scale m GUT ~ 1.2 x 10 16 GeV. This allows to re- 
duce the number of free parameters to few, and calculate 
the rest by performing renormalization group evolution. 
Such approach gives a much more convincing picture of 
the low energy spectrum than just assessing the running 
masses of the superparticles. We have also demonstrated 
the influence of the mixing among squarks and quarks on 
the results. The evident weak point of our approach was 
the lack of knowledge of the trilinear coupling constants, 
which forced us to use additional source of constraints 
(neutrino oscillations, 0^2/3 decay). 

In this paper we study the Majorana neutrino transi- 
tion magnetic moments generated in the i?-parity violat- 
ing modified minimal supergravity model. This model, 
based on the Ref. [l9| , takes into account the full depen- 
dence of the renormalization group equations (RGE) on 
the RpV couplings and their soft SUSY breaking ver- 
sions. All these couplings have been incorporated in 
the unification scheme, so, within this assumption, the 
whole model can be calculated numerically. No other in- 
put, in particular experimental data, is needed. We con- 
centrate on the contributions coming from the charged 
Higgs bosons, charged higgsinos, charged gauge bosons, 
charginos, leptons, and sleptons. The neutral particles, 
like the neutral components of the Higgs bosons, neu- 
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tralinos and sneutrinos will contribute to the neutrino 
mass, but will not contribute to the magnetic moment, 
and as such will not be discussed here. 

The paper is organized as follows. After the intro- 
duction, in Sect. II the minimal supersymmetric stan- 
dard model with i?-parity violation and supersymmc- 
try broken by the gravitational interactions is presented. 
Sect. Ill is devoted to the technical problems of han- 
dling the free parameters and calculating the coupling 
constants in the Higgs sector during the electroweak sym- 
metry breaking. The Majorana neutrino transition mag- 
netic moments, and the interesting for us contributions 
to this quantity are discussed in Sect. IV. At the end 
numerical results are presented. 



II. ^-PARITY VIOLATING SUGRA 

The minimal supersymmetric standard model (MSSM) 
[2ll [22[ is the minimal extension of the standard model 
(SM) of particles and interactions, which introduces su- 
persymmetry (SUSY). In short, for each particle there 
is a new superpartner (fermionic for bosons and bosonic 
for fermions) so that the number of fermionic and bosonic 
degrees of freedom in the MSSM are equal. The one ad- 
ditional multiplet, not present in the SM, is a second 
Higgs doublet needed to cancel the gauge anomalies and 
to give masses to both the up and down components of 
the SU(2) doublets. This version of the model preserves 
the so-called i?-parity, which is a multiplicative quantum 
number defined as R = (— i) 2S + 3 ( B -- L ) 5 5 being the spin 
of the particle, B the baryon, and L the lepton number. 
One easily sees, that R = +1 for ordinary standard model 
particles, while R = — 1 for the supersymmetric partners. 
Preserving R not only gives the baryon and lepton num- 
ber conservations but also forbids decays of SUSY parti- 
cles into SM particles, leaving the lightest SUSY particle 
stable. 

The problem of the lepton and baryon number con- 
servation is still disputable. On one hand, our observa- 
tions suggest, that both B and L are conserved quantum 
numbers. Especially strong limits on B violation come 
from the non-observation of the proton decay. On the 
other hand, in the SM not only the total lepton num- 
ber but also generation lepton numbers L e , L^, and L T 
are conserved, and this rule has been proven wrong after 
the discovery of neutrino oscillations. From the theo- 
retical point of view, there is no underlying symmetry 
or mechanism, which supports conservation of B and L, 



which might suggest that it is only an accidental (not 
exact) symmetry present in the low-energy regime. It 
is generally expected, that the higher-energy extensions 
of the SM will not preserve at least the lepton number 
(the proton stability must still hold). Such models based 
on supersymmetry introduce new interactions, or in fact 
do not rule out certain terms, that should be present in 
the superpotential, which violate i?-parity |23l427| | . They 
have a much richer phenomenology and predict lots of 
new phenomena, like the lepton number violating neutri- 
noless double beta decay, RpV loop-corrected neutrino 
mass and magnetic moment, and others. 

We work within the slightly modified model described 
in details in Ref. [l9j . It is defined by the following R- 
parity conserving and i?-parity violating parts of the su- 
perpotential 



where 



RpC 
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n(A*E k )ijLfLjEk + (A D k)ijLfQj l 'D kx 



Cub 

le xyz (Aui) jk U?D]Dt - e abK l L«H b u . (3) 



Here Y's are the 3x3 trilinear Yukawa-like couplings, 
fi the bilinear Higgs coupling, and A and k % represent 
the magnitudes of the i?-parity violating trilinear and 
bilinear terms. L and Q denote the SU(2) left-handed 
doublets, while E, U and D are the right-handed lepton, 
up-quark and down-quark SU(2) singlets, respectively. 
Hd and H u mean two Higgs doublets. We have intro- 
duced color indices x, y, z = 1,2,3, generation indices 
k = 1,2,3 = 6, /1, t and the SU{2) gauge indices 
a,b = 1,2. Finally, e ab and e xyz with e i2 = £123 = 1 
are the totally antisymmetric tensor densities. 

The supersymmetry is not observed in the regime of 
energies accessible to our experiments. Therefore it must 
be broken at some point. A convenient method to take 
this fact into account is to introduce explicit terms, which 
break supersymmetry in a soft way, ie., they do not suffer 
from ultraviolet divergencies. We add them in the form 
of a scalar Lagrangian [19j, 



- C = m 2 Hd h\h d + m\h\h u + lHm 2 L )l + k\m 2 L%Hd )h d + h d (m 2 HdL Jk + q\m 2 Q )q + e(m 2 E )e^ + ^(m^)^ + u(m^)^ 

+ i (Mx&B + MzWiW* + M 3 g~Jg a + h.c.) (4) 

+ [(A E )ijkh d e.j + : A/, !}.,(!, + (Au)ijqih u Uj - Bh d h u 

+ (A E k) i jlil j e k + (A D k)ijliqjdk + {Aui)j k Uidjdk - D,kh u + h.c], 
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where the lower case letter denotes the scalar part of the 
respective super-field. Mi are the gaugino masses, and A 
(B, Di) are the soft supersymmetry breaking equivalents 
of the trilinear (bilinear) couplings from the superpoten- 
tial. 



III. FREE PARAMETERS AND THE 
LOW-ENERGY SPECTRUM 

The main problem of most supersymmetric models is 
their large number of free parameters. Therefore addi- 
tional constraints are usually imposed, which introduce 
certain relations among them, effectively reducing their 
number. It is a well known fact, that the RGE equations 
for the gauge coupling constants lead in the MSSM model 
to their unification at energy jtigut ~ 1-2 x 10 16 GeV. 
This feature is absent in the SM, although suggested by 
the extrapolation of the LEP1 data. It may therefore 
seem natural to assume also certain type of unification 
of other parameters. Let us introduce at mcuT- 

• the common mass of all the scalars mo, 

• the common mass of all the fermions TOx/2, 

• the common proportionality factor Aq for the soft 
SUSY breaking couplings 



squark mass eigenstates as 



W,D,E 



— AnY 



1 U,D,Ei 



A ui>D i tE i=A Q Y u , D , E (i = 1,2,3). (5) 

This scheme of unification follows the mSUGRA assump- 
tions with the exception that we do not assume the total 
universality of the A couplings, but vary them keeping 
them proportional to the Yukawa Y coupling constants 
at the unification scale. As a free parameter remains 
also the ratio of the Higgs vacuum expectation values, 
tan/3 = v u /vd, and the sign of the /i coupling constant, 
sgn(/i). The sixth free parameter in our considerations is 
the initial value of the A's at mz, A . We leave this value 
free and investigate its influence on the results. The idea 
is that we want to have A's non-zero and contributing 
from the beginning to the RGE running of other param- 
eters. In this way, Ao controls the amount of i?-parity 
violation in the low energy regime. Notice, that A's will 
not unify at mcuTi althou gh t heir RGE running is al- 
most flat (cf. Fig. 1 in Ref. |20j|). 

The procedure of obtaining the low energy spectrum of 
the model involves few renormalization group runnings 
between the low (mz) and high energy scale (mcuT), 
and between mz and the scale at which the electroweak 
symmetry breaking occurs, which is defined by the top 
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At this scale the radiative corrections to the Higgs po- 
tential coming from the squarks are minimal. We start 
by fixing the Yukawa couplings at mz using the quark 
and lepton mass matrices Mjj d e 



M v 
M D 
M E 



vS Dr Y^ Dl , 



(7) 



where S matrices perform diagonalization so that one 
obtains eigenstates in the mass representation. The 
Yukawa-like RpV couplings at mz are all set to 



A 1, 



(8) 



1 being a 3 x 3 unit matrix. Next, we evolve the gauge 
and Yukawa couplings from mz to the unification scale. 
During this running we turn off the Higgs sector, /1 — 
Ki = B = Di = 0, and set the sneutrino vev's Vi = 0. 
These parameters will be calculated later. 

At mcuT we impose the GUT conditions as described 
above, ie., we set the massess of all the scalars to mo, of 
all the fermions to m 1 / 2 , and the soft breaking couplings 
to AqY. The Yukawa couplings Y's are left unchanged 
in order to reproduce correctly the masses of leptons and 
quarks when evolved down to mz- So are the A's, which 
are free parameters here. We evolve all the running pa- 
rameters, coupling constants and masses, down, and find 
the best scale for minimizing the scalar potential Eq. © 
and at this scale calculate \i = sgn(/^) yj/IF and B using 



m 



1/4 



m 2 Hd — to 2 ^ tan 2 j3 M\ 



tan 2 /3 - 1 



B 



sin 2/3 



(m 2 Hd - m 2 Hu + 2\n\ 2 ). 



(9) 
(10) 



Running all the parameters down q m ; n — > mz non-zero /j, 
and B generate non-zero Ki and Di . These are used as the 
starting values to the second long run. Next, we perform 
once again the RGE running up to the wqut scale, but 
now with all the Higgs parameters contributing to the 
running. We impose GUT conditions, go down, calculate 
the new q m - m scale and once again evaluate the corrected 
Higgs parameters, this time together with the sneutrino 
vev's, according to 



tan 2 (3-1 



Z 1 * * 

m UH d — + K if 1 — 



I'd 



I'd 



'H, 



+ «i 



V u 



M 2 



tar//n--^, (11) 
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B = 



sin 2/3 



+ (ml Hd +K*^-D 



I'd 



(12) 



ui[(m^)n + |ki| 2 + D'] + w 2 [(m|) 2 i + Ki«£] + u 3 [(m|) 3 i + kikI] = -[m^ dii + + 

«i[(mi)ia + /t 2 «a] + «a[(mi)aa + l«2| 2 + D '] + t> 3 [(m|) 32 + k 2 k 3 ] = -[™?H d L 2 + P*Ka]vd + D 2 v u , (13) 

u i[( m i)i3 + k 3 k*] + v 2 [(m 2 L ) 23 + k 3 k 2 ] + "3[K) 3 3 + |k 3 | 2 +-D'] = -[m 2 Hd L 3 + M*«3]«d + D 3 u„, 

where D' = Mf 22 ^ + (g 2 + flf) S ir^( u2 - v 2 , - vj). The tadpole equations (O can be easily solved and we get for 
the sneutrino vev's 



deW 



1,2,3, (14) 



where 



(ml) 11 + |ki| 2 + D' (m 2 L ) 21 + Ki k* 2 (m 2 L ) 31 + Ki k* 3 

W = ( (m 2 L ) 12 + K2 kI (m 2 L ) 22 + \k 2 \ 2 + D' (m 2 L ) 32 + k 2 k* 3 |. (15) 

(m 2 ) 13 + K 3 K* (m 2 L ) 23 + K 3 n* (m 2 ) 33 + |k 3 | 2 + D' 

and Wi can be obtained from W by replacing the i-th column with 



-[m 2 ?dL . 2 + H*K 2 \Vd + X> 2 V n )• (16) 
[ m2 H d L a + M* K 3]«d + £>3«« 



The newly obtained v^s are incorporated in the scheme 
in the following way: 

v u = sm/3^v 2 - (vf + v% + uf), (17) 
u d = cos(3^Jv 2 - (v 2 +« 2 + t;3), (18) 

f 2 = (246 GeV) 2 , which preserves the definition of the 
tan P and w 2 = w 2 + v \ + v\ + v 2 + v 2 . Notice, that in this 
model sneutrino vev's contribute to the mass of the Z 
and W bosons. The equations (|ll[) - (fT3")) are solved sub- 
sequently until convergence and self-consistency of the 
results is obtained. In practice three iterations turned 
out to be enough. After this procedure we add also the 
dominant radiative corrections |28j |. and get back to the 
mz scale to obtain the mass spectrum of the model. 



IV. NEUTRINO MAGNETIC MOMENTS 

It is well known that in the RpV supersymmetric mod- 
els only one neutrino gains mass after the diagonalization 
of the neutrino-neutralino mass matrix. The remaining 
contributions come from the sneutrino vev's at the tree 
level and from the one-loop processes. The main one-loop 
mechanism involves decomposing a Majorana neutrino- 
neutrino vertex into a particle-sparticle loop, which basi- 
cally is either the quark-squark, or the (charged) lepton- 
slepton loop. These contributions will be proportional to 
some functions of the masses of the particles and (A^>) 2 



or (Ac) 2 . Detailed discussions can be found, e.g., in 
Refs. 0411 El- 

Neutrinos are neutral particles and as such cannot in- 
teract with photons. However, this interaction is possi- 
ble through the one-loop mechanism, in which charged 
particles appear in the loop. Therefore, an effective neu- 
trino magnetic moment may be generated. Due to the 
CPT theorem, Majorana neutrinos cannot have diagonal 
magnetic moments, which act between the same flavours 
of neutrinos, but they may have the off-diagonal tran- 
sition magnetic moments [J]. To be more precise, the 
transition Majorana magnetic moment acts between vu, 
and Uj L chiral components of Majorana neutrinos, as- 
suming gauge theory with only left-handed neutrinos. As 
a consequence it violates the total lepton number by two 
units (AX = 2), and can be discussed provided that R- 
parity violation occurs. The effective Hamiltonian takes 
the form 

#cff = l^L<J aP ^ L F Pa + h.c. (19) 

The contributing Feynman diagrams are presented on 
Fig. [TJ We will evaluate the amplitudes of these pro- 
cesses in the interaction basis, in which the vertex cou- 
pling constants are known directly from the mixing ma- 
trices. However, one has to take into account, that the 
propagators of particles travelling inside the loop must 
be written for their mass eigenstates. This can be done 
using the mixed-propagator formalism [29j . Assume the 
following relation between the mass (J) and interaction 



(a) 



(b) 



(c) 



(d) 



(c) 



FIG. 1: Feynman diagrams with two mass insertions inside the loop contributing to the Majorana neutrino transition magnetic 
moments. 



(a) eigenstates of some particle, 



integral, so it may be written in a general form as 



(20) 



Then, the amplitude in the weak basis can be written in 
terms of the Green function and expanded in the mass 
basis 



(21) 



so the flavour changing amplitude may be written as 

A(a^8)=\ \ 11 (22) 
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A more detailed derivation and discussion of this result 
can be found in Sect. 4.2 of Ref. [29j . 

The matrix V is obtained during the diagonalization 
procedure of the mixing matrices and its columns are the 
eigenvectors of the mixing matrix. The magnetic moment 
will consist of the product of numerical constants, cou- 
pling constants, and the function 1 describing the loop 



2m e fj, B , (23) 



where N c is the color index (=3 for quarks and squarks, 
=1 otherwise), are the trilinear coupling constants, 
Xi y2 are the mass insertions residing inside the loop, Q 
is the charge of the particle the photon is attached to 
(in units of the elementary charge), m e is the electron 
mass (2m e rj 10 -3 GeV), and hb the Bohr magneton. 
The sum runs over all mass eigenstates which form every 
mixed particle present in the loop. 

The masses and couplings in our approach are calcu- 
lated from the boundary conditions at the mz and toqut 
scales (see Sect. III). The loop integrals X can be de- 
rived analytically using standard techniques of integra- 
tion in the Minkowski space. Denoting by /i,/2, • • • and 
61,62, ■• • the fermions and bosons which appear in the 
loop, and by m,t } b their respective masses, the loop inte- 
grals for diagrams (a)-(e) depicted on Fig. [T] read 



l{a) = m h m h ^(/l, hM, h) +-7 r 4(/i,/2,&i, h) + ™, fl mf 2 ^(/a, h, h, h) + ^(h, fi, 61, 6 2 ), 

J< 6 > = m h -F 3 (/i,6 1 ,6 2 ,6 3 ), 

1^ = 4 J- 4 (6 1 ,/ 1 ,/ 2 ,/ 3 ), 

l(d) = m h m h ( m h + m h ) Fa (/a > A > h , h ) + (2m fl + 3m /2 + m h ) T 4 (f 2 , fi , / 3 , h ) , 

x[e) = m fA m h m h + m f2 m f3 + m / 3 m /i)- ;r 3(/i,/2,/3,&i) + (5m/! +2m/ 2 +2m/ 3 ) Ti{fuf2,hM) 

+ m h ( m /i 171 f 2 + m h 171 h + m f 3 m fi)J 7 3(f3jij2,b 1 ) + (5m/ 3 + 2m/ 2 + 2m fl ) J" 4 (/ 3 , h , f 2 , 61 ) 

where the functions are given by 

d 4 k 



(167T 2 )J" 3 (a,6,c,d) = (16tt 2 

m l log(m 2 ) 



1 

(2tt) 4 (k 2 - m 2 a ) 2 (k 2 - m 2 )(k 2 - m 2 c )(k 2 - to 2 ) 
r^-r - l) - 1 



mz —mz 



ml log(m 2 ) 



(24) 
(25) 
(26) 
(27) 

(28) 



(m 2 a - m 2 )(m 2 a - m 2 c )(m 2 a - to 2 ) 
m 2 c log(TO 2 ) 



(m 2 -to 2 )(to 2 - m 2 c )(m 2 



ml) 



(167r 2 )J-4(a, 6, c, d) 



(to 2 — to 2 )(to 2 — m 2 )(m 2 — to 2 ,) (mz, — Tn1)(m 2 d — m 2 )(m 2 d — to 2 ) ' 

2 r°° _d^k k 2 

J-00 W) 1 (fc 2 - m 2 a ) 2 (k 2 - m 2 )(k 2 - m 2 )(k 2 - to 2 ) 



(29) 



6 



mf t \og(m%) ( — 

a o\ a; \ m ^_ 7 



mt —m. 



m\ log(m^) 



™b log(mg) 



(mg - m 2 a )(m 2 b - m 2 c ){m 2 - m 2 d ) 



mj \og(m 2 d ) 



(m 2 - m 2 a )(m 2 c - m 2 )(m 2 - m 2 ) (m 2 d - m 2 a ){m 2 d - m 2 )(m 2 - m 2 ) 



(30) 



We have kept denotations from Ref. [l3| . Please notice 
also that a typo has been corrected in Eq. (|25|) . and new 
terms appeared in Eqs. (IM1) and (|2"5)) . 



V. RESULTS 

In order to contribute to the neutrino magnetic mo- 
ment, the loop must contain charged particles. The 
quarks-squarks and leptons-sleptons have already been 
discussed elsewhere (see Sect. I). The discussion on 
the bilinear insertions on the neutrino lines, possible 
due to the neutrino-neutralino mixing, can be found 
in Refs. [3, G3]- I n this work we have focused on the 
charged Higgs-sleptons and the leptons-charginos inter- 
action eigenstates. 

Altogether we have identified 45 diagrams of the forms 
given on Fig. [T] containing the particles in question. The 
way of obtaining them is a straightforward but rather 
tedious task, and we have used a computer program to 
match all possible vertices and internal lines to the given 
patterns. The trilinear vertices can be constructed di- 
rectly from the superpotential ((T|), and the possible bi- 
linear mass insertions are defined by the mixing matrices, 
as given in Ref. fl9j . 

A few conditions must have been met for the result 
to be accepted. First of all, during the RGE evolution 
some of the particles may get negative mass parame- 
ters squared. Such tachyonic behaviour resulted in im- 
mediate rejecting of the given input parameters. Also, 
finite values for the Yukawa couplings, which tend to 
'explode' for too low or too high values of the tan/3, 
were required. When analyzing the results, we have re- 
jected points, for which the masses of the superparti- 
cles where too low. The mass limits we have used where 
taken from the Particle Data Group report [3Cj , and they 
read: m^o > 46 GeV, m^o > 62 GeV, m^o > 100 GeV, 
m^o > life GeV, m-± > #4 GeV, m-± > 9 3 4 GeV, m s > 
107 4 GeV, m A > 94 GeV, m? > 82 GeV, m§ > 379 GeV, 
nig > 308 GeV. This is by far the conservative choice 
of constraints. The new results from the LHC projects 
push the quoted above limits to higher values, closer to 
the 1 TeV range. We do not use them here for various 
reasons. Firstly, most of the results are still marked as 
preliminary, and more careful data analysis is being per- 
formed right now. Also, the LHC data are interpreted 
within the simplest SUSY models, with many additional 
assumptions, like the absence of i?-parity violation, and 
equality of the gluino and squark masses. Theses as- 
sumptions are not true in the model used here. 



A. Constraining the model's parameter space with 
the magnetic moments 

There are different assessments and constraints put on 
the neutrino magnetic moment. Analyzing the impact 
of the solar neutrino data on the neutrino spin-flavour- 
precession mechanism, the authors of Ref. [3l| have found 
a model-dependent upper value for the magnetic moment 
to be fewx 10 _12 /is. On the other hand, direct searches 
in the MUNU experiment resulted in the upper limit to 
be 9x 10 -11 /ib [32|]. In this paper we adopt a conservative 
limit and reject points, for which \x v > lO -10 //^. A more 
strict approach would result in a substantially narrower 
parameter space of the model. 

The results are presented on Figs.[2][9l for \x > 0, A = 
10~ 4 , A) = 200 - 1000 GeV in steps of 100, m . 1/2 = 
200 - 1000 GeV in steps of 20, and tan/3 = 5 - 40 in 
steps of 5. Each point represents a valid set of input 
parameters, which results in a physically accepted low 
energy spectrum. Each figure represents a fixed tan/9, 
and panels read in rows from left to right correspond to 
A taking values 200 GeV, 300 GeV, . . . , 1000 GeV. 

In general, there is no global regularity in these results. 
Certain regions are excluded due to unaccepted values of 
the masses. The excluded stripes, visible in several of 
the panels, come from too high values of ijl v . Starting 
from tan/3 = 15, a parabolic-like shape, whose position 
depends on the Aq parameter, appears. It starts to be 
visible for A — 400 GeV and crosses the to — TO1/2 
plane with increasing Aq. [We are not sure about the 
origin of this 'stability line'. As a side remark we add, 
that it appeared also in our analysis of the masses of 
the lightest Higgs bosons h° in the context of the CDF- 
DO discovery [33|]. Constructing similar plots, but using 
the condition, that 120 GeV < m h o < 140 GeV, we have 
observed very similar parabolas.] We notice also, that for 
tan /3 = 40 the points look totally random, which shows 
that the model breaks down in this region due to too high 
values of the Yukawa couplings. 



B. Ao dependence 

Another interesting problem is how the initial value 
of the RpV couplings at mz, represented by a common 
parameter A , influences the results. We have chosen to 
check this relation for a few specific sets of parameters. 
For this, the Snowmass SUSY benchmark points [34| were 
used. The dependence of the resulting contribution to the 
transition magnetic moments have been calculated for 
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FIG. 3: Like Fig. [2] but for tan j3 = 10. 
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FIG. 4: Like Pig. [2] but for tan /3 = 15. 
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FIG. 5: Like Fig. [2] but for tan j3 = 20. 
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FIG. 6: Like Fig. [2] but for tun/3 = 25. 
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FIG. 7: Like Fig. [2] but for tan j3 = 30. 
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the SUGRA benchmark points SPSla, SPSlb, SPS2, and 
SPS3, arc shown on Figs. [TUHT51 In this calculation we 
have not discriminated the results which were exceeding 
the upper limit for 

The existing limits on A's, obtained withtin SUSY 
models, oscillate around 10 -2 — 10 -5 , depending on the 
method used [ll|, HU . We check here the range between 
fewxlO" 2 and 10~ 10 . We notice that the RGE running 
decreases the values of A's for higher energies [20j, thus 
the parameter Ao sets the upper limit on them. 

We first notice, that below roughly Ao = 10~ 5 (in some 
cases even more) the results stabilize and do not change 
with decreasing Ao, which makes them almost indistin- 
guishable with the situation when all RpV couplings are 
set to zero. On the other hand, the calculations broke 
down for Ao > few x 10~ 2 . This leaves a rather nar- 
row region of Ao, between 1 and 3 orders of magnitude 
and only modestly exceeding 10~ 2 , for which the RpV 
couplings play a role. 

We see also that for small Ao the general alignment 
of the /x's resembles hierarchical structure, while for high 
Ao this hierarchy vanishes, and the results show no reg- 
ular pattern. Too high values of Ao can boost them 
to values close to 1Q~Vb, which are excluded by ex- 
periments like MUNU [32J (recall upper limit from the 
MUNU ~ 10~ 10 /is). Curiously, for small (or zero) val- 
ues of Ao, all the transition magnetic moments tend to be 
of the order of 10~ 10 fi B (SPSla, SPSlb) or lower (SPS2, 



SPS3). 

Let us also compare the newly computed contributions 
with similar contributions coming from simplest loops 
with no mass insertions (trilinear RpV couplings only) 
and loops with bilinear insertions on the external neu- 
trino lines (neutrino-neutralino mixing). We do it for 
two points, for which earlier calculations have been pre- 
sented in Refs. [H], [H| - the unification with low val- 
ues: Ao = 100 GeV, mo = m 1 / 2 = 150, and high values: 
A = 500 GeV, m = m 1/2 = 1000. In both of these 
cases tan j3 = 19 and fi > 0. 

The results are shown in Tab. Q] for the low point 
and Tab. |TT] for the high unification point. Below the 
numerical values, results for pure trilinear loops [161 ] and 
bilinear contributions are given as ranges, obtained 
for different cases (assumptions of the normal or inverted 
hierarchy of neutrino masses, data from the neutrinolcss 
double beta decay searches). 

For the LOW point the smallest contributions discussed 
in this paper are at least 5 orders of magnitude greater 
than any other calculated so far. This clearly shows, 
that the loops with heaviest particles (like higgsinos and 
charginos) tend to dominate the overall contribution to 
the magnetic moments. This result is not totally un- 
expected, but it was not fully clear, if the higher order 
process (amplitude proportional to four instead of two 
coupling constants) will not suppress the effect coming 
directly from the masses of the heavy particles. That 
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FIG. 9: Like Fig. [2] but for tan/3 = 40. Please notice, that for such high value the model breaks down and the results are no 
longer reliable. 
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VI. SUMMARY 



was the case when the neutrino-neutralino mixing placed 
two mass insertions on the external neutrino lines of the 
diagrams. Here, however, the mass dependent loop func- 
tions I are capable to boost the amplitudes of the process 
to very high values. All this holds also for the second, 
HIGH, set of input parameters, although there the differ- 
ence between smallest value of fi Vl2 and the maximum 
trilinear contribution is 'only' two orders of magnitude. 



The full contribution to the Majorana neutrino transi- 
tion magnetic moment consists of three main parts. The 
first one is represented by loop diagrams with two tri- 
linear RpV couplings, containing either a quark-squark 
or lepton-slepton pair. The second takes into account 
possible neutrino-neutralino mixing, which occurs on the 
external neutrino lines in the form of bilinear mass inser- 
tions. The third one, discussed in this paper, allows for 
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trino for SPS3 SUGRA point. 

the mass insertions to appear inside the loop, which may 
contain a number of different particles, i.e., charged Higgs 
bosons and the corresponding higgsinos, leptons and slep- 
tons, charged gauge bosons, and charginos. A proper dis- 
cussion of the latter contribution is possible only within 
a consistent i?-parity violating model, in which all the 
mixing between different mass eigenstates is taken into 
account. Also, new terms in the RGE equations, propor- 
tional to the RpV couplings appear. However, the impact 
of these terms has already been studied [2(| and their 
presence changed the low energy mass spectrum of the 
model by a factor of 1/5 at most (20%). This indicates, 
that the crucial point is the phenomenon of mixing be- 
tween different mass eigenstates, so that the amplitudes 
of the processes must be expanded in the physical bases 
and summed over respective eigenstates. 

We have found that the presently discussed contribu- 
tions are dominant over the remaining two. This allowed 
us to find the acceptable parameter space of the model, 
using the condition that no magnetic moment may exceed 



10~ fin- We have also checked, how different initial val- 
ues of the RpV couplings, represented by Ao, change the 
results. Finally, a comparison with previous numerical 
studies has been given. 

TABLE I: Comparison of the magnitudes of different contri- 
butions to the Majorana neutrino transition magnetic mo- 
ments for the input parameters set LOW (see text). Here sub- 
scripts 1, 2, 3 = e, /i, r. 
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TABLE II: Comparison of the magnitudes of different con- 
tributions to the Majorana neutrino transition magnetic mo- 
ments for the input parameters set high (see text). Here 
subscripts 1, 2, 3 = e, /i, r. 
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